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ABSTRACT 



Establishing the stellar masses, and hence specific star-formation rates of submillimetre galaxies is crucial for deter- 
mining the role of such objects in the cosmic history of galaxy/star formation. However, there is as yet no consensus 
over the typical stellar masses of submillimetre galaxies, as illustrated by the widely differing results reported from 
recent optical-infrared studies of submillimetre galaxies with spectroscopic redshifts z ~ 2-3. Specifically, even for the 

Osame set of submillimetre galaxies, the reported average stellar masses have ranged over an order of magnitude, from 
~ 5 x 10 10 Mq to ~ 5 x 10 11 M©. Here we study how different methods of analysis can lead to such widely varying 
• results. We find that, contrary to recent claims in the literature, potential contamination of IRAC 3-8 fj,m photome- 

try from hot dust associated with an active nucleus is not the origin of the published discrepancies in derived stellar 
masses. Instead, we expose in detail how inferred stellar mass depends on assumptions made in the photometric fitting, 
and quantify the individual and cumulative effects of different choices of initial mass function, different "brands" of 
Q ■ evolutionary synthesis models, and different forms of assumed star- formation history. We review current observational 

S_t 1 evidence for and against these alternatives as well as clues from the hydrodynamical simulations, and conclude that, for 

^ | the most justifiable choices of these model inputs, the average stellar mass of luminous (Ssso > 5mJy) submillimetre 

galaxies is ~ 2 x 10 1 M© to within a factor ~ 2. We also check and confirm that this number is perfectly reasonable 
in the light of the latest measurements of the dynamical masses of these objects (~ 2 - 6 x 10 11 M from CO (1-0) 
observations), and the evolving stellar mass function of the overall galaxy population. Galaxy stellar masses of this order 
imply that the average specific star-formation rate of submillimetre galaxies is comparable to that of other star-forming 
galaxies at z > 2, at 2-3 Gyr -1 . This supports the view that, while rare outliers may be found at any stellar mass, 
00. most submillimetre galaxies simply form the top end of the 'main-sequence' of star-forming galaxies at these redshifts. 

Conversely, this argues strongly against the viewpoint, frequently simply asserted in the literature, that submillimetre 
galaxies are extreme pathological objects, of little relevance in the cosmic history of star-formation. 
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1. Introduction 



X 



2010bt iRicciardelli et all 120101: lHavward et al] |2011allbL 
20121 ). Unless the uncertainty in stellar masses can be re- 



The properties of galaxies selected at submillimetre wave- ducedj it is in essence impossible to tell whether SMGs 
lengths (SMGs) have yet to be fully understood. In par- have specific s tar-formation rates (sSFR = SFR/M*) corn- 



ed ■ ticular, there is still no real consensus over their stel- parab i e to those of other star-forming galaxy populations 

lar masses. A robust measurement of the stellar masses at z > 2j or whether they are extreme outliers from the 

of SMGs is key to addressing a number of important ' ma in-sequen ce' of star-forming gala xies j Daddi et alJ SOul 

issues, including the evolution of the relationship be- iNoeske et al]l2007l IGonzalez et alJl201Ch . 
tween galaxy and black -hole masses ([Borvs et al.l I2005t 

lHainline etld]|2011|, [mi| thereafter), and the contribution Most previous studies of the stellar masses of SMGs 

of SMGs to the overall h i story of stellar mass assembly have con cluded in favour of values in the range ~ 10 11 - 

(|Michalowski et al.| BoTfli, |M10| hereafter). Perhaps most 10 12 M Q (iBorvs et al.ll2005b iDve et al.ll2008f ). However, two 

importantly, a proper understanding of the stellar masses recent stud ies of ~ 70 SMGs wit h spectroscopic redshifts 

of SMGs is required to distinguish whether their massive ~ 2-3 from lChapman et al.l (|2005l) have reached apparently 



dust-enshrouded star-formation rate s are fuelled primarily very different conclusions, highlighting t he impact of alter- 
by mergers or by cold gas infal l (ISw nbank et al.l l2008t native methods of analysis. Specifically, iMichatowski et al.l 



iDave et al.| |2010fc IGonzalez etTal] |201lfc |Narayanan et al] (|2010al fbh found stellar masses consistent with previous 



studies , with a median value M» = 3.5 x 10 n M R (see also 



Scottish Universities Physics Alliance lHatsukade et al.l l2010t llkarashi et al.l 1 2 1 it ISantini et al.l 
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l2010t iTamura et aT1l2010t) . whereas lHainline et all (|201lh 

have reported values, based on the same SMGs, which are 
systematically lower by a factor of ^ 6 (m edian M* = 5- 
7 x 10 10 M o ; see also IWardlow et al.ll2011fl . 

The primary objective of this paper is to explore the 
origins of this discrepancy, and to attempt to resolve it. 
More generally, however, we have taken this opportunity to 
expose in detail how the derived stellar masses (and hence 
sSFRs) of SMGs depend on different assumed forms of star- 
formation history, and different "brands" of evolutionary 
synthesis models. The use of a sample of SMGs with spec- 
troscopic rcdshifts is especially useful for this purpose, since 
it is well known that, if only photometric rcdshifts arc avail- 
able, the errors in derived stellar masses can be dominated 
by the effect of uncertainties in z p hot- 

The structure of this paper is as follows. In Scction[5]wc 
delineate the various factors and choices of analysis meth- 
ods which influence the stellar mass determinations and 
quantify their impact on the final derived values. We also 
review recent evidence that informs the preferred choice 
of, for example, initial mass function (IMF) and evolution- 
ary synthesis model, and conclude that for most reasonable 
choices of such model inputs the derived stellar masses in 
fact only differ by a factor ~ 2. In Section [3] we investi- 
gate the consistency of our re-derived stellar masses with 
the latest measurements of the dynamical masses of these 
objects, the evolving stellar mass function of the overall 
galaxy population, and the predictions of theoretical sim- 
ulations. Wc discuss the implications of our results on the 
sSFR of SMGs in Section 0] and close with a summary of 
our conclusions in Section [5] We use a cosmological model 
with H Q = TOkms^Mpc- 1 , fl A = 0.7 and O m = 0.3. 

2. Determining the stellar masses of SMGs 

2.1. Initial mass function 

The choice of a particular IMF introduces a systematic un- 
certainty of a factor of ~ 2 in the determination of the stel- 
lar m asses (M*) and star-formation rates (SFRs: lErb et al.l 
|2006|) . Sp ecifica l ly, the differen ces b etwe e n the IMFs as - 
sumed by IMIOI (ISaloeterl Il955h and Ell dGhabrierl l200l 
accounts for a factor of 1.8 in the difference between the 
quoted stellar mass estimates 

Both top-heavy and standard IMFs have been claimed 

*h et al.l 
This 

confusing situation is further exacerbated by the fact 
that different studies of local ellipticals (arguably the 
likely present-day descendants of SMGs) have yielded 
a wide range of p r eferred IMFs i n cluding bottom-light 
(Ivan Dokkuml 120081). iKroupal (1200 lh (iGerhard et all 1200 it 
ICappellari e t al.l I2006T). Salpeter dGrillo et all I2009D. and 



dulii Lup-iieavy cinu sLanuciiu iivir s nave ueen cj 
to reproduce the number counts of SMGs dBaugl 
120051: iFontanot etahl l2007t lHavward et ail 1201 lbl ). 



eve n bottom-heavy (Ivan Dokkum fc Conrovl20iaEHll: sec 



also lThomson fc Charvll201ll ). 

Hence, since there is as yet no clear evidence that the 
IMF of SMGs is systematically differ ent from that of other 
galaxies, in this paper we assume the lChabrierl {2003) IMF 
(unless stated otherwise) to which all appropriate derived 
qua ntities are conv erted if necessary (values derived using 
the lSalpeter!ll955l IMF are divided by 1.8). This IMF has 
been derived for the Milky Way and represents an interme- 
diate choice between the top-heavy and the bottom-heavy 
IMFs. We also note that this "Chabrier IMF" is essentially 




Fig. 1. Spectral energy distribution (flux density vs. ob- 
served wavelength) of SMMJ123606.85+621021.4 at the 
spectroscopic redshift z = 2.509 fitted using the 

double-compone nt star formation history method of 

ICirasuolo et al.l (|2010l . IClOf ) . Black circles denote the data, 
red and blue curves are the SEDs of the old and young 
components, respectively, and the green curve is the sum 
of these two components. The older component dominates 
the near-IR light and the derived stellar mass, whereas the 
young component delivers almost all rest-frame ultraviolet 
flux density. 



identical to th e "Can onical IMF" recently summarised by 
IWeidner et~aTI (|201lD . While we argue that this IMF is the 
natural choice on the basis of current evidence, it must still 
be accepted that we must live with a fundamental uncer- 
tainty of ~ x2 in stellar masses until the uncertainty in 
the IMF in SMGs is resolved. However, we also note that 
the derived sSFR is unaffected by this IMF uncertainty, 
making it a particularly useful quantity for comparing the 
properties of different types of star-forming galaxies. 

2.2. Evolutionary models and star-formation histories 

It is well known that the choice of the evolutionary syn- 
thesis model used to fit the optical-infrared photometry 
can affect the derived stellar masses of any type of galaxy. 
The G RASIL spectral energy distribution (SEP) tem- 
pl ates (Silva et all 119981 : llglesias-Paramo et al. 20071) used 
by MIC make use of the Padova stellar tracks (lAlongi et al 



1 19931: iBressan et al.lll993l:lBertelli et alJl994H F agotto et al 
ll994aUbl ld: lGirardi et al.lll996ll2000r i. The input metallicity 
for GRASIL was calculated based on the chemical evolution 
model in which individual stellar populations are allowed 
to have metallicity between 0.2-5 Zq. The mean (median) 
me tallic ity d erived for SMGs is 8 Zq (0.6 Z e ). 

IHUI used iBruzual fc Charlotl (|2003ft models (which are 
also based on the Padoya libra ry) and i n add ition explored 
the use of the Maraston| (|2005l) models. iHllI assumed solar 
metallicity in order to derive a typical mass-t o-ligh t ratio. 

The stellar masses derived by [mi] using 
IBruzual fc Charlotl (|2003l) models are o n average ~ 50% 
higher than those calculated using the IMarastonl (|2005l) 
models. This average difference is at least in part a con- 
sequence of the higher contribution of thermal ly pulsing 
asymp totic giant branch (TP-AGB) stars in the IMarastonl 
(|2005l ) models, which reduces the derived mass-to-light 
ratio at rest-frame near-infrared wavelengths at stellar 
population ages ~ 1 Gyr (a point discussed further below). 

An apparently separate (but, it transpires, in fact con- 
nected) cause of differing stellar mass estimates is the cho- 
sen paramcterisation of the star-formation histories (SFHs) 
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Table 1. Comparison of stellar masses and specific star formation rates of SMGs calculated using various methods, 
stellar population models and star formation histories. 



Fit" 


SSP" 


SFH C 


Mean Af« 


Median M, 


Mean SSFR 


Median SSFR 


Age Y d 


Ageo e 


Fraco 1 








(log Mq) 


(logM Q ) 


(Gyr" 1 ) 


(Gyr" 1 ) 


(Gyr) 


(Gyr) 




M10 


Padova 


cont+ burst 


11.32 ±0.05 


11 31+ uuv 

— 0.06 


2.86 ±0.43 


1.74±»-i5 


0.050 


1.856 


0.84 


Hll 


BC03 


lburst /const 


10.87 ±0.06 


10.88_ ; 07 


9.71 ± 1.99 


4.60l ;g 


0.820 






Hll 


M05 


lburst /const 


10.71 ±0.06 


10 74+ ' 08 
iu. f4_ .05 


12.96 ± 2.17 


a or + 1.56 
D.OO_2.09 


0.860 






CIO 


BC03 


double 


11.44 ±0.08 


11.471°;- 


10.87 ±8.68 


i oc+0.21 
1 - oo -0.31 


0.191 


1.332 


0.87 


CIO 


BC03 


tau 


11.16 ±0.10 




17.36 ±11.98 


2 56+ 1 - 00 


1.178 






CIO 


BC03 


lburst 


11.06 ±0.10 


ii n2 +0 06 


14.88 ± 8.69 


3.07l i 6 7 


0.246 






CIO 


M05 


double 


11.73 ±0.06 


1 1 79+O.O3 


1.25 ±0.19 


o.67tg:?; 


0.075 


1.777 


0.66 


CIO 


M05 


tau 


10.46 ±0.13 


10.78+n 7 7 


167.34 ± 47.42 


r i fi +6.38 


0.767 






CIO 


M05 


lburst 


11.24 ±0.07 


11 22+ 012 


8.43 ± 4.57 


1 qn +0 - 49 
i.yu_ 43 


0.198 







Notes. The IChabrieri (|2003f ) IMF is assu med. Errors are standard deviations for the means and bootstrap 68% ranges for 
medians. (a) SED fitting method presented in lMld (|MichaIowski et alJl2010ah . lHlll l|Hainline et al.ll201ll "). or lClOl ijCirasuolo et all 
I2010T 1. (6) Single stellar population models being either Padova tracks, BC03 (|Bruzual fc CharlodfeOOSf ) or M05 (|Marastonll2005ri . 
' c ' Assumed star formation history being either a continous SFH with a burst (cont±burst), an average of the single burst and 
constant (lburst /const), double burst (double), single burst (lburst), or exponentially declining (tau). ^ Mean age of the younger 
stellar component. ' e ' Mean age of the older stellar component. Mean fractional contribution to the total stellar mass of the 
older component. 



of the galaxies in question. As noted bv lHlll due to their 
high redshift, and often rather red spectral-energy distribu- 
tions, it is very difficult to r ecover robust SFHs for SMGs 
from photometric data (cf. iDve et al.l [20081 12010ft . This 
leads inevitably to difficulty in constraining galaxy ages and 
mass-to-light ratios. 

In an attempt to account for this uncertainty, IH11I ex- 
plored the use of two alternative single-component SFHs, 
adopting cither an inst antan eous starburst or a continuous 
star-formation history. iHUl then averaged the two result- 
ing estimates to pro duce the adopted stellar mass for each 
galaxy. In contrast, IMIOI assumed a two-component SFH 
consisting of a smooth Schmidt-type law (SFR oc M gas ) 
with a recent starburst superimposed commencing 50 Myr 
before the observed galaxy SED was emitted, and with a 
mass in the range 0-9% of the infalling gas mass. 

The assumption of a multi-component SFH generally 
naturally leads to higher mass-to-light ratios and, in turn, 
higher stellar masses than the use of a single-component 
model. This is because, while the starburst component can 
account for the ultraviolet (UV) emission (and the infrared 
emission from the UV- heated dust), the second (older) com- 
ponent is then free to contribute more stars with higher 
mass-to-light ratios (see Fig. Q] and also Schael et al. in 
preparation). By contrast, the use of a single instantaneous 
starburst model limits the age of the entire stellar popula- 
tion to the young age of the starburst required to match the 
UV emission (and thus the true stellar masses are inevitably 
undcr-cstimatcd), while in the continuous star- formation 
model the current SFR is set by the UV emission and the 
total age is then limited in order not to overshoot the op- 
tical and the near-IR part of the spectrum (assuming the 
galaxy has always formed stars at the same high rate). 

In order to further clarify the impact of the choice 
of stellar population synthesis models and star-formation 
histories on the derive d stellar masses we ha ve revisited 
the SED fitting of the IChapman etall ([20051 ) SMG sam- 
ple. Specifically , we have u sed the method presented in 
ICirasuolo et al.1 (120071 I2010D. based on the Hyp erz pack- 
age utilising both lBruzual fc Chariot! (|2003ft and IMarastonl 



(|2005D models and assuming three different forms of SFH: 
an instantaneous burst of star-formation, an exponentially 
declining star- formation rate (so-called tau models), and a 
double-component model composed of two instantaneous 
bursts with different ages and (independent) dust attenu- 
ations (an example is shown on Fig. [T] and all estimates 
are listed in Tab. I A. l]) . In the double-component model the 
age of the young component was varied between 50 Myr and 
1.5 Gyr, and the old component was allowed to contribute 
0-100% of the near-IR emission while its age was varied 
over the range 1-6 Gyr. Solar metallicity was assumed. 
These models were run using the Chabrier IMF. 

The results of this modelling are illustrated in Fig. [5] 
and summarised in Table [IJ where the mean and median 
values of M* and sSFR (as well as ages and contributions 
of the older population) for the SMGs are tabulated for the 
full range of different stellar populatio n syn t hesis models 
and star- formation histories explored in IMIOI . IHUl and the 
present work (it is not possible to distinguish the models 
based on the goodness of fit). 

What is immediately evident from this table is that 
the results are strongly dependent on the assumed pa- 
rameterisation of the star-formation history (see also 
lBussmannetlII[20T2T) . as is the impact of assuming differ- 
ent population s ynthesis m o dels. I n particular, it is note- 
worthy that the IMarastonl (|2005ft models only result in 
significantly lower stellar masses (and hence substa ntially 
higher sSFRs) than the iBruzual fc Charlotl (|2003ft mod- 
els when a single-component exponentially-declining tau- 
model SFH is assumed. We investigated the origin of this 
effect, and found that, when allowed only single-component 
tau models, the mean preferred age of the IMarastonl (|2005ft 
model fits is ~ 0.8 Gyr (see Tab. [T]), which corresponds 
closely to the time at which the en hanced TP-AGB star 
contribution in the IMarastonl (|2005h models has maximum 
impact. This result provides fu rther evidence th at the com- 
bination of tau SFHs and the [MarastorJ ([20051 ) models is 
not appropriate for SMGs. 

Reassuringly, however, there is no significant difference 
between the masses derived with these two alternative mod- 
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.0 10.5 11.0 11.5 12.0 
log Stellar mass (ivy - MichatowskMO 



10.5 11.0 11.5 12.0 
log Stellar mass (ivy - Hainlinel 1 



Fig. 2. Comparison of the stell ar masses of SMGs as der ived bv iMlOl. iHllI and in t he present work using the stellar 
population synthesis models of iBruzual k Chariot! (|2003, toy ) and iMarastonl (|2005l bottom). The different symbols 
indicate the stellar masses derived by IMlOl (re d circles), Hill (blue stars), and our new estimates based on the SED 
modelling method of lCirasuolo et ail <|2010l IClOh for three alternative forms of star formation history: double component 
model (red squares), exponentially declining tau mo del (light blue triangles) and single burst (green plu ses). A ll the derived 
masses have been converted to the lChabrierl (|2003l) IMF. Small symbols denote z < .5 SM Gs. The lMlOl stellar masses 
are consistent with our new estimates for the double-component models, whereas the IHllI stellar masses are consistent 
with single component models, and arc on average a factor of ~ 2.5 smaller than the double-component estimates (see 
Tab. □]). 



els for a single burst SFH. This is because, in order to re- 
produce the star-forming nature of SMGs, the single burst 
fitting typically results in much younger ages of ~ 0.2 Gyr, 
corresponding to the time when the TP-AGB contribution 
is small and does not have a significant influence on the 
derived stellar masses. 

As expected (see above discussion of the lMlOl modelling 
approach), the stellar masses obtained using the double 
componen t SFH are systematically larger. Especially, the 
IMarastonl (|2005f ) models with the double SFH yield stellar 
masses higher by ~ 1 dex than these derived assuming the 
tau SFHs with these models. 

In summary, the derived stellar masses depend primarily 
on the assumed parameterisation of the SFH, and the dif- 
ferent evolutionary synthesis models only yield significantly 



different results if the single-component exponentially- 
declining SFH is adopted (due primarily to the extreme be- 
haviour of the Maraston 2005 models with this adopted pa- 
rameterization of SFH). Since, as we argue below, a model 
of a rapidly starforming galaxy in which the SFR is assumed 
to have been even higher at all previous epochs is obviously 
physically unrealistic for SMGs, the preferred choice of evo- 
lutionary synthesis model is actually of less importance in 
deriving the stellar masses of SMGs than might have been 
supposed. Nevertheless, in the next section, we do, for com- 
pleten ess, review the a rgu ments for and again s t the choice 
of the IMarastonl (|2005ft or IBruzual k Charlotl (|2003l ) mod- 
els. 

Of more significance is the fact that, as can be seen 
from Table [TJ our new two-component models using the 
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iBruzual fc Charlotl (l2003h m odels yield a value very simi- 
lar to that derived by Mid using a different form of two- 
component model, and a different evolutionary synthe- 
sis model (albeit based on the same stellar tracks as the 
IBruzual fe Charlotll2003l models'). Thus, consistent average 
stellar masses are produced by models which allow recent 
and past SFR to be decoupled, regardless of whether the 
earlier star-formation is parametrised as a single burst or a 
more extended period of star formation. This consistency 
is both reassuring and helpful given that, as we argue fur- 
ther below, a SFH history allowing potentially decoupled 
recent and past star-formation activity is rather obviously 
appropriate for objects selected on the basis of recent very 
violent star formation within an already metal enriched en- 
vironment. 



2.3. Choice of evolutionary synthesis model 

Recent years have seen some vigorous debate over whether 
the stronger TP-AGB cont ribution in t roduc ed into evolu- 
tionary synthesis models bv lMarastonl (|2005l ) is justified by 
the data. 

Some worker s have argued for the superiority of the 
IMarastonl (120051) models based on derived galaxy stel- 
lar masses (Ivan der Wei et abl 2006aB ). on derived photo- 



metric redshi fts dMaraston et al 



-II2006D. on average broad- 
band SEDs (|Maraston et all 120061: ICimatti et all 120081) 



on o ptical and near-infrared colours ([MacArthur et all 
[20T1. and bv c ompar ison with se mi-analytical models 
(|Henriques et all 1201 ll ). Moreover, IZibetti et all (|2009h 
have claimed that stellar mass determinations based on 
i- and ii-band luminosities are more consistent with each 
other when long-lived TP-AGB stars are introduced. 

Ho wever lEminian et al.l (|2008| ) showed that IMarastonl 
(120051 ) models predic t exce ssively red optical colours, and 



recently iKriek et al.l (|2010h have found that the average 
spectral energy distribution of post-starburst galaxies at 
0-7 z ^ 2.0 does not appear to display the excess near- 
inf rared flux predic ted from the TP-AGB contribution in 
the lMarastonl (I2005T) models. I ndeed, they demonstrate that 
the IBruzual fc Charlotl (|2003j ) model provides a much bet- 
ter description of the the data, even at the age when the 
TP-A GB contribution should be near maximum. In ad- 
dition, IConrov fe Gunnl fl201dl their Fig. 2 and 16) have 
found that the IMarastonl (20051) models predict excessively 
red colours for post-starburst gala xies and for star clu sters 
in the Magellanic Clouds. Finally. IZibetti et al.l <|2012h did 
not find any strong carbon feature of TP-AGB stars in the 
near-IR spectra of post-starburst galaxies. 

The work presented here cannot be used to definitively 
resolve this issue, but we note that the extreme outlying 
nature of the tau-model IMarastonl (120051 ) results shown in 
Table Q] arguably provides some further evidence against the 
plausibility and appropriateness of the strong TP-AGB con- 
tribution adopte d in these models for tau SFHs. Moreover, 
we note that the IMarastonl (|2005| ) models with the double 
SFH yield the highest median stellar masses for SMGs. 

We also note tha t, confining t he comparison to the use o f 
the lChabrierl (|2003l ) IMF and the lBruzual fe Charlotl (120031) 
models, the median stellar mass of SMGs derived bv lMlOl 
(M* = 2 x 10 11 M©) is in fa ct only a factor of ~ 2.5 higher 
than that reported bv lHllI (M* = 7.6 x 10 10 M o ). 



2.4. Choice of star-formation history 

The key issue for stellar mass determination, therefore, ap- 
pears to be the choice of how to parameterize the SFH of 
SMGs. Reassuringly, however, the results do not appear to 
be very sensitive to precisely how this is parameterized, pro- 
vided the SFH allows past and "present" star-formation ac- 
tivity to be, at least in principle, disconnected. This can be 
seen from the fact that the different two-component models 
based on the Padova tracks listed in Table Q] yield average 
stellar masses in the range M* = 2-3 x 10 11 M Q . 

It only therefore remains to consider if these two- 
component models should be preferred to single-component 
burst or tau models. In fact, at least for SMGs, this is a 
relatively straightforward argument, since both the single- 
burst and tau models both enforce the condition that the 
star-formation rates in SMGs were always higher at epochs 
earlier than that at which the object has been selected 
due to its extreme star-formation rate. Indeed this prob- 
lem lias_ajrea^y_b^endiscussed by several authors includ- 
ing [MarastorreTjiLl (|2010l ). They concluded that, for high- 
redshift galaxies, the assumption of an exponentially de- 
creasing SFH (as well as constant SFH to some extent) 
leads to unrcalistically low ages and under-estimated stel- 
lar masses, because t he SEP fit is complete ly dominated 
by young stars. In fact iMaraston et all (|2010l) found an ex- 
ponentially increasing SFH could provide a much better 
description o f data . Such SFHs are conceptually similar to 
those used in[M10, with the SFR being highest close to the 
epoch at which the observed galaxy SED was emitted, and 
indeed there is growing evidence in favour of rising SFHs 
in the wi der population of g alaxies studied at hig h red- 
shift (e.g. lFinlator et alJl20ll1: iPapovich et al.ll201lli The 



lack of eit her tw o-compon ent models or rising S FHs in the 
analysis of iHllI (or indeed IWardlow et al.ll201lh can there- 
fore explain all of the rema ining discr epan cy in the average 
stellar masses reported by IMIOI and IHllI since, as can be 
seen from Table [TJ the use of a double component SFH 
results in stellar masses of SMGs being ~ 3 times larger 
than that deduced using the single burst model and ~ 2 
times larger than that derived using the tau models. This 
is also illustrated in Figure [2 where the masses derived 
assuming single-component SFHs (light blue triangles and 
green plusses) are systematically lower than those derived 
assuming the more realistic double-component SFHs (red 
squares). This conclusion is consistent with the findings of 
Schael et al. (in preparation^ who found that th e stellar 
masses of the SHADES SMGs dCoppin et al.ll2006h derived 
using a two-component SFH are on average a factor of 
~ 2 — 3 higher than when derived using a single-component 
SFH. 



2.5. Other model-fitting issues 
2.5.1. Wavelength range 

The wavelength range used in the SED model fitting could 
potentially affect the resulting derived stellar masses. For 
completeness we therefore tested whether t he diff eren ce in 
the SMG stellar mass estimates derived by IMIOI an d IHllI 
is not, at least in part, due to the fact that iMlfl fitted 
the entir e mu lti-frequency dataset up to radio wavelengths, 
whereas IHllI confined attention to observed A < 8.0 /im. 



See also lSchaell (|2009h . 
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Table 2. Stellar, gas and dynamical masses of SMGs (in units of 10 11 M©) 



SMG 


M*(Mio) a 


A^*(cio) 6 


M*(H11) C 






. ^-highCO d 
iW dyn 


M CO(1-0) e 


SMMJ123549.44+621536.8 


9 

1. 00 


+2.65 
-1.33 


U.00 


+0.66 
-0.33 


2.09+"'f 7 

— U.4 ( 


9 H7 


± 


1.73 J 


U.90 


n A A Q 

± U.44 


2.30 ± 0.40 1 


SMMJ123618.33+621550.5 


2.78 


+2.77 
-1.39 


4.27 


+4.25 
-2.13 


0.76±°o1 


0.35 


± 


0.11 s 


2.18 


± 0.18 s 


4.36 ± 0.36' 


SMMJ123634. 51+621241.0 


2.54 


+2.53 
-1.27 


0.71 


+0.70 
-0.35 


0.87+1;? 


0.25 


± 


0.08 s 


1.12 


± 0.25 s 


2.24 + 0.51' 


SMMJ123707.21+621408.1 


4.95 


+4.93 
-2.47 


1.12 


+ 1.12 
-0.56 


-i 09+0.32 

1 - OZ -0.27 


1.69 


± 


1A1 S 


2.78 


± 1.05 s 


2.20 ± 0.60 7 


SMMJ131201. 17+424208.1 


1.27 


+ 1.27 
-0.63 


6.76 


+6.73 
-3.37 


1 02+ 33 

1 - uz -0.25 


1.34 


± 


0.50 s 


9.50 


± 2.38 s 


19.01 ± 4.76 ! 


SMMJ163650.43+405734.5 


2.06 


+2.05 
-1.03 


3.80 


+3.78 
-1.90 


1 ' 4O -0.32 


2.56 


± 


2.14 / 


3.50 


± 1.79 s 


4.10±0.80 / 


SMMJ163658. 19+410523.8 


3.27 


+3.26 
-1.63 


8.13 


+8.09 
-4.05 


1 oq+0-33 


2.95 


± 


2.46 / 


1.40 


± 0.88 s 


5.80 ± 0.40 / 


SMMJ163706.51+405313.8 


4.52 


+4.49 
-2.25 


10.23 


+ 10.18 
-5.10 


1 62+ 024 

1 - oz -0.21 


0.24 


± 


0.07 h 


3 


.40^ 


6.80* 


SMMJ221735. 15+001537.2 


1.30 


+ 1.30 
-0.65 


1.91 


+ 1.90 
-0.95 


>S1+ 06 

u - dl -0.06 


0.30 


± 


0.07' 1 


2 


.80 h 


5.60 1 



Notes. (a) lMichalowski et all (|2010al ) with lChabr icr ( 2003) IMF. {b) From this work using the method of lClOl for lBruzual fc Charlotl 
((20031 ) models and double SFH. (c) lHainline et al] |201ll ) for iBruzual fc Charlotl (pOuD models. (d) Dynamical mass from high 
exitation CO line. (e) Dynamical mass from CO(1-0) line. (/) llvison et all (|2uTll ). (s) lEngel et al.l (|2010l ). (h) ICreve et all ((20051 ). 
^ Obtained by multiplying the M^^ lCO by a factor of 2, the average correction from high excitation lines to CO(1-0) derived for 
four SMGs with CO(1-0) data. Mass errors for lMlOl andlClOl are simply a factor of 2. 



To te st this we simply repeated the analysis performed by 
iMlfl but this time restricting the model fitting to pho- 
tometry at optical, near-infrared and IRAC wavelengths. 
We obtained almost identical results with a median stel- 
lar mass of ~ 1.7 x 10 n M Q instead of ~ 1.9 x lO n M , 
and hence conclude that the longer-wavelength data have 
no significant influence on the derived stellar masses. This 
is reassuring and is easily explained since the emission at 
far-infrared wavelengths depends only on the properties of 
the dust and of the younger stellar population that is heat- 
ing the dust, and should be essentially unaffected by the 
older stars which dominate the total stellar masses of SMGs 
((Mini) . 



2.5.2. AGN contamination 

SMGs are believed to suffer only a minor (~ 10-30%) 
contamination of their bol ometrie luminos i ties by ac- 
tive galactic nuclei (AGNs: lAlexander et al.l 120051 | 2008f 



Menendez-Delmestre et al.l 2007. 2009; Valiante et al.l2007 



Pope et al.ll2008t lHainline et alJl2^fMlMurphv et al.l l2009; 
Watabe et all l2009L ~ iMlOh . Here we consider whether the 
presence of AGNs in SMGs can significantly influence the 
derivation of their typical stellar masses. 

iHUl reported that, for ~ 50% of SMGs, non-stellar 
(AGN) power-law features contribute < 10% in the rest- 
frame ii-band, whereas, for ~ 10% of SMGs, the contri- 
bution is > 50%. However, even for those objects which 
do appear to have a significant AGN contribution at ob- 
served IRAC wavelengths, it is not obvious that the stellar 
masses based on our full multi-wavelength fitting will be 
significantly distorted, especially given the typically large 
(~ 20%) photometric errors generally assigned to IRAC 
data in the SED fitting p rocess . In fact we find that the 
stellar masses reported by IHUl are consistent with our es- 
timates for single component models (Fig. [5] and Tab. [TJ, 
without making any corrections for possible AGN contribu- 
tions. We therefore do not find any evidence t o supp ort the 
idea that the lower stellar masses reported bv lHUl are due 
to the subtraction of an AGN component. Rather, as de- 
scribed abo ve, it seems clear that the lower stellar masses 
deduced by IHUl arise, primarily, from the adoption of a 
single-component SFH. 



We also note that the presence of any power-law features 
in the IRAC dat a sho uld not have strongly affected the 
mass estimates of IMIOL simply because the SED fits were 
not forced to reproduce the power-law features and, in fact, 
the best-fitting stellar models tend to under-cstimate the 
longer wavelength IRAC data when the p ower- law compo- 
nent appears to be present (see Fig. Al of lMlOh . Moreover, 
we do not fin d any correlation between the power-law frac- 
tio n give n bv lHlll and the stellar masses or M*/Lk derived 
bv lMlOt which would be expected if the power-law features 
were biasing the lMlOl stellar mass estimates high. 

2.5.3. Lensing 

If SM Gs are lensed then the stellar masses derived by 
iMldl and IHUl should be corrected down. At l east some of 



1998: Downes & Solomon 20031: 


iDunloD et al.l 2004: 


Kneib et al.l 12004 iMotohara et al. 


120051 Wilson et al.l 


2008: Berciano Alba et al. 20101: 


Knudsen et al. 20101: 


Neerello et al.l 2010t Swinbank et a 


[. I2010bt IVieira et al.l 


2010t iTamura et al.l 20101 jHezaveh & Holder] 2011: 


Ikarashi et alJ 20111). However. Lima et al.l ((20101 Fie. 2) 



found that the number counts of SMGs with Ss5o < 20 
mjy do not require any l e nsing correction. Moreover, the 
modelling of iPaciea et al.l (|2009l . Fig. 6) showed that only 
< 10% of SMGs with Ssso < 15 mJy at z < 3 are magnified 
by a factor of > 2. Hence, the lensing correction for SMGs 
in the 850 /im flux-density range considered here is likely 
to be small on average. 

3. Consistency tests for stellar masses 

3.1. Dynamical masses 

The estimates of the dynamical and H2 gas masses from 
CO lines detected for SMGs provide a consistency check on 
the stellar mass estimates. Specifically, assuming that the 
extent of the CO-emitting gas is the same as the extent 
of the stellar component, the sum of the gas and stellar 
masses h as to be less t han t he dynamical mass. Using this 
approach lEngel et al.1 (|2010t) derived a mean ratio (AI* + 
M?as)/- Mdvn — 2.1 for s even SMGs usin g the stellar masses 
of lMlOl (corrected to the lChabriei1l2003l IMF) and concluded 
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t hat the stellar masses derived bv lMldl are inconsistent with 
their gas and dynamical masses. As explained below, we do 
not find that this is the case. 

There is, of course, some uncertainty and contro- 
versy over the calculation of dynamical a nd molecular 
gas masses from CO li ne observations (e.g. iGenzel et al.l 
120031 : llvison et aDl201ll ). Perhaps of most relevance f or the 
present discussion is the fact that lEngel et"al"1 ([2010D used 
high excitation lines (CO J=7-6, 6-5, 4-3 and 3-2), which 
trace denser and hotter gas, very likely confined to the more 
central regions of a galaxy. Therefore, the derived dynam- 
ical masses could easily refer to regions smaller than the 
extent of the stellar component, making it difficult to com- 
pare them directly with the ste llar masses. Indeed , based 
on low excitation CO(l-O) lines llvison et all (|2011l) found 
dynamical masses higher by a factor ^ 2.4- 4.1 for two out 
of three SMGs, for which lEngel et all (|2010t) claimed a sig- 
nificant discrepancy between stellar and dynamical mass 
(H DF 76 and N2 850.2; the third one was not observed 
by llvison et all ). Moreover, recently iRiechers et al.l (| 2 1 If) 
fo und CO(l-O) sizes ~ 60-90% larger than that reported 
bv lEngel etldl (|2010l ) for two additional SMGs. 

Interestingly, the mean CO half-light radius obtained 
bv lEngel et atflj2010l) of 2.4 kpc is indeed 15-40% smaller 
than the typical half-light radius of the near-infrared emis- 
sion dominated by th e stellar component (2.8-3.4 kpc; 
Swinbank et~aTI l2010at iMosleh et all 120111: iTargett et all 



1 10 
Dynamical mass (10 11 ivy 



Fig. 3. Stellar mass (M*) versus dynamical mass (Md yn ) 
estimated from CO(1-0) (upper panel), and stellar+i^ gas 
mass (M* + A/ ga s) versus Md yn (lower panel) for the 9 
SMGs for which the necessary results have been compiled 
in Table [3J The d ata points derived from the stellar mass 
estimates of iHllI are indicated by blu e stars, those based 
on the stellar mass estimates of iMlOl ar e ind icated by red 
circles, while the new double-component ICld values are in- 
dicated by red squares. In principle, only the area below 
the dotted line represents a consistent set of galaxy prop- 
erties with A/g as + M* < Mdyn, but because of the typical 
scatter (indicated by the error bars in the lower panel) and 
unknown inclinations, ~ 2 — 3 out of the 9 individual SMGs 
are in fact expected to exceed this apparent threshold (see 
Section I3~I] for details ) . We see that, contrary to the con- 
clusion of lEngel et al.l (l2010h. t here is no evidence that the 
stellar masses derived bv lMlOl using the Chabrier IMF are 
in conflict with the dynamical mass estimates, even allow- 
ing for estimates of the additional gas mass. Indeed, the 
average values of these ratios favour the IMlOl stellar mass 
estimates as the most consistent with the available data 
(see Section l3~Tj) . 



201 it) . Mor e over, one out of three SMGs for which 



Engel et al.1 (|2010l ) claimed inconsistency with the M10 



stellar mass (N2 850.2) was observed bv ISwinbank et al 
(|2010al ) who measured the stellar extent to b e ~ 3 times 
larger than the extent of the CO(7-6) found bv lEngel et al.l 
(I2010D . This reinforces the likelihood that the dynami- 
cal masses of the problematic SMGs have been under- 
estimated, as the stellar component is more extended than 
the high-excitation gas. 

This argument is supported by hydrod y namical 
simulations of SMG s (|Naravanan et al.1 120091 l2010bl: 
lHavward etHI l2011allbT ). where only ~ 30% of the stel- 
lar mass is distributed ove r the inner r adius of 2.5 kpc, 
roughly the size over which lEngel et al.l (|2010l ) calculated 
dynamical masses. This correction to the dynamical masses 
would be eno ugh to remove the inconsistency with the stel- 
lar masses of Mil even though some simulations predict 
that the usual conversions result in the dynamical masses 
overp redicted by a factor of ~ 1.5-1.9 ([Narayanan et al.l 
I2009L Fig. 10 and Eq. 6). 

In an attempt to overcome these systematics, and facil- 
itate a fairer comparison of dynamical and stellar masses, 
we focus here on the newly available CO (1-0) dynamical 
masses which, as described above, appear to better trace 
the full extent of the SMGs, yielding somewhat larger sizes 
and masses. Stellar, gas and dynamical masses are avail- 
able for 9 of the gal axies in our sample , including two CO 
detections made by iGreve et al.l (|2005l ) . These dateJEI are 
presented in Tabled 

In the final column of this table we give the CO (1- 
0) dynamical mass as our best estimate of the dynamical 
mass of each galaxy. Where available (4 SMGs) this num- 
ber has been taken directly from the CO (1-0) observations 



2 Note that lEngel et al.l (l2010f) used the stellar masses re- 
ported initially in Hainlind pOOa ) whereas we quote the cor- 
rected values from IH 111 . 



Michalowski et al.: The stellar masses and specific star- formation rates of submillimetre galaxies 



performed bv llvison et al.l (|201lh . while for the remaining 
objects this number has been inferred from the higher-order 
CO-line dynamical masses by multiplying them by the av- 
erage correction factor for the af orementioned 4 SMG s for 
which the dynam ical masses of llvison et al.l (|2011[ ) and 
lEngel et al.l (fcoiOt) can be directly compared. 

In Fig. [3] we then plot, for these 9 SMGs, M* and also 
Mgas+M* as a function of Md yn , with the dynamical masses 
taken from the final column in Tabled and the stellar mass 
estimates f rom jHlll iMldl and the present work (double- 
component [Clflmodel) indicated by the different symbols. 

From Fig. [3] it can be seen that, within the errors, there 
is no evidence that even the higher stellar mass determina- 
tions for any object exceed the estimated dynamical mass. 
At first sight it might seem that there could be more of a 
problem (for all mass determinations) for M* + M gas , but 
this is not really the case, both because of uncertainty in 
M gas , and because of the potential effects of varying in- 
clination angle on Md yn - Specifically, while the unknown 
inclination (i) of gas discs should not influence the average 
(M gas + Af * ) /Mdyn ratio, if the real inclination of an indi- 
vidual galaxy is lower than the assumed average, then its 
derived Md yn is inevitably too low, which can make it ap- 
pear to be inconsistent with the (orientation-independent) 
derived stellar mass. To quantify this effect we consid- 
ered 100,000 galaxies with real (M gas + M»)/M dyn = 1 
(or 0.5) with inclinations distributed uniformly in three di- 
mensions, i.e. with a mean (sin 2 i) = 2/3. We found that 
the adoption of this mean value for all galaxies results 
in ~ 18% (9%) of them appearing to exceed the allowed 
(M gas + M»)/M dy „ = 1 by a factor of > 2 and ~ 12% (7%) 
by a factor of > 3. Hence, in a sample of ~ 10 galaxies 
with random unknown inclinations, we would expect 1-2 
of them to significantly exceed the allowed ratio, simply 
because of the unknown inclination of any molecular disc 
in an individual object. 

In summary, one should expect, given the effects of ran- 
dom inclination and uncertainty in the measurement of M* , 
Mgas and Md yn , that a few galaxies should lie above the 
one-to-one relation indicated by the dotted line in Fig. [31 
and it is the mean (or median) value of the mass ratios 
which need to be considered when reviewing the viability 
of the stellar masses. Working from the numbers given in 
Table ^ for the ratio M*/Mdv„ the mean (median) val- 
ues are 0.36 (0.24) for the Ell] stellar masses, 0.80 (0.64) 
for iMlOl and 0.94 (0.78) for the new doub le-component 
IClOl values. Thus, adoption of the lower iHlH stellar masses 
would imply that only ~ 1/3 of the dynamical mass of a 
typic al S MG c an be accounted for by M*. In contrast, the 
IMlOl and IClOl results would imply that M* for a typical 
SMG is ~ 3/4 x M dy „. If the values of M gas in Table H 
are taken at face value then this latter scenario would ap- 
pear to be more consistent with the data, as for the ratio 
(M* + M gas ) / Mdvn the mean (median) values are 0.71 (0.50) 
for the |H11| stellar masses, 1.15 (1.07)for [MH and 1.09 
(1.19) for the new double-component IClOl values. This then 
prov ides further evidence that the st ellar mass estimates of 
iHlli are a little too low, the new lClOl do uble- component val- 
ues are perhaps a little high, and the IMlOl results appear 
to provide the best description of the data. 

We conclude that the latest estimates of the dynamical 
masses of SMGs are in very good accord with the stellar 
masses of SMGs estimated byMl3 with the lChabrierl (|200l 
IMF. 



3.2. Number density of massive galaxies 

Another consistency check on the plausibility of the stellar- 
mass estimates is to assess how the number density of SMGs 
compares with the number density of galaxies in general 
in the same implied stellar-mass range, at both lower and 
comparable rcdshifts. 



iDrorv et all (|2004l) reported the cumulative number 
density of galaxies n(M* > 1.1 x 1O U M ) = 3.16- 
0.64 x 10~ 4 Mpc~ 3 at redshifts z = 0.5 - 1.1 (their Fig. 6 
and Table 3, after conversion of their mass limit to the 
IChabrierll200a IMF) . Restricting the analysis to the S MGs 
at z > 1.4, and using the comoving volumes given bvlMlOl 
we obtain n SMG {M* > 1.1 x 1O U M ) = 7.0 x 10" 6 Mpc- 3 
for the M* v alues of IMlOl and 3.8 x 10 ~ 6 M pc~ 3 for the 
M* values of IHlli . Hence, even for the IMlOl stellar mass 
estimates the number density of SMGs is ~ 30 times lower 
than that of massive galaxies at lower redshifts and there- 
fore the number of massive SMGs is clearly not inconsistent 
with the number density of massive galaxies at lower red- 
shifts, even allowing for a relatively modest and plausible 
duty cycle (see below) . In addition, from the mass compar- 
isons discussed in the previous subsection, further conver- 
sion of the observed molecular gas reservoir into stars, even 
if ex tremely efficient, can at most increase the measured 
MIOl stellar masses by ~ 30 — 50% (or, equivalently, the 
Hill masses by ~ 100%), and so this basic conclusion is not 



seriously affected by potential further star formation. 

The redshift range 1.4 < z < 3.6 corresponds to a 
time interval of 2.7 Gyr. Given the typical gas reservoirs 
of M gas — 1 x 10 11 M Q listed in Table[2]it can be seen that 
the typical gas exhaustion timescale is ~ 200 Myr (assum- 
ing a typical SFR of ~ 500 M yr -1 ). If one further assumes 
that a typical SMG is observed mid-way through its final 
violent star-forming peak, then a reasonable estimate of the 
sub millimetre luminous lifetime of an SMG is 400 Myr (see 
also lBouche et al1l2005l) . This implies a duty-cycle number 
density correction factor of ~ 6 — 7 for the redshift interval 
in question. If one further factors an incom pleteness correc- 
tion of a factor of ~ 3.5 as derived bv lMlOl . then the implied 
true number density of potential SMGs in this era is close to, 
but still less than the best estimate of M* > 1.1 x 10 11 M Q 
galaxies at z ~ 1.1. Hence, the conclusion that the stellar 
masses of SMGs are typically M* ~ 2 X 10 11 M appears to 
be perfectly consistent with the number density of massive 
galaxies at z ~ 1. 

A similar conclusion can be reached comparin g to the 
nu mber density reported by iPozzetti et al.l (|2007l Fig. 11) 
and lKaiisawa et all (|2009l Fig. 16): n(M„ > 6xlO lo M ) ~ 
10~ 3 -2 x 10~ 4 Mpc~ 3 at z = 0-2, whereas the correspond- 
ing values for SMGs are tismg(M. > 6 x 10 10 M Q ) = 
8.7 x 10~ 6 and 5.6 x 10~ 6 Mpc" 3 . 

Final ly, we integra t ed the mass function at z = 3-3.5 
given by (Caputi et all (|201lD for M* > 10 11 M galaxies, 
obtaining 7 x 10~ 6 Mpc~ 3 . This rather narrow redshift in- 
terval corresponds to an elapsed time of only 350 Myr, and 
so, following the above arguments, the duty cycle correc- 
tion should be ~ unity. It is therefore perhaps unsurprising 
that this number density is rather similar to the number 
density of SMGs. 
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Fig. 4. Submillimetre flux densities as a function of the 
stellar masses of SMGs. Solid lines represent the evo- 
lutionar y tracks from the merger hyd ro dynamic simu- 
lations (|Naravanan et all I20091 l2010al lbl; lHavward et al.l 
l2011al fbT) enclosed in the yellow shaded area, the big square 
represents the median stellar mass of 2 x 10 10 Mq for 
SMGs with 5Wj_> 5 mJy in the semi-a n alytical models 



dBaugh et al.l 120051: 



1201 It I Almeida et al 



Swinbank et al.l 120081 : iGonzalez et al.l 
12011 ). whereas t he big circle corre- 



sponds to the median value derived by iDave et all (|2010h 
from the cosmological hydrodynamic simulations. Other 
symbols represent the o bservational data for SMGs with 
stellar masses from lMlOl (red circle s), [Hill (blue stars) and 
our new estimates based on ICld and double-component 
SFH (red squares). Small symbols indicate SMGs at z < 
0.5. According to the hydrod ynam ic simulations, many of 
the stellar masses derived by iHlli arc too low to result in 
measur able su bmillimetre flux, while the stellar masses de- 
rived bv lMld are consistent with the hydrodynamic simula- 
tio n prediction s (with in the observational errors). Neither 
the EH IMKH or [C10| stellar masses appear to agree, on av- 
erage, with the median stellar mass predicted by the semi- 
analytic models. 



3.3. Numerical simulations of SMGs 

Figure [4] shows submillimetre flux density versus stellar 
mass for a set of idealized hydrodynamic simulations of 
merging and isolated disk galaxies, where 3-D radiative 
transfer has been performed to calculate the submillime- 
tre flux^JTlie_j3hnulat ions are_similar to those described 
Narayanan et al.l (|2009L l2010afbh : see lHavward etaD 



for a description of the differences. Note that 
the plot simply shows tracks for different simulations. No 
weighting by abundance has been applied, so only the range 
spanned by the simulations, not typical values, can be read 
from the plot. 

We note that in order to reach 5*850 > 5mJy, 
these simulations require stellar mass to exceed M* ~ 
1 x 10 11 M Q . As can be seen from Fig. |H the predicted 
stellar mass range is i n goo d agreement with the SMG stel- 
lar masses derived bv lMlOl We note that, in the cold infall 
model of SMGs, comparab ly high stellar masses are also 
preferred ([Dave et aLlboiOL their Fig. 2). 



Gonzale z et al.l ( 2011 



ha nd, th e semi-analytic al mod - 
(|2005h . ISwinbank et all (I2008T) 
Fig. 3) and lAlmeida et al.l (f201lh 
predict significantly lower stellar masses for SMGs, with 
a median Af* ~ 2 x 10 10 M (indicated by the square on 
Fig. 2]). Such a value is lower even t han the vast majority of 
th e stellar masses derive d bv IHlli confirming the findings 
of ISwinbank et al.l (|2008| ) that the stellar masses of SMGs 



are generally systematically under-predicted by this class 
of semi-analytic model. 

4. Specific Star-Formation Rates 

One important step towards locating SMGs within the gen- 
eral framework of galaxy/star formation at high redshift is 
to establish whether, on average, they simply represent the 
most massive extension of the normal star-forming galaxy 
populatio n in the young ac tive Universe at z > 2 (as pre- 
dicted by IDave et ail 120101 in cosmological hydrodynamic 
simulations), or whether they are extreme outliers in terms 
of specific star-formation rate (sSFR). 

To explore this we compare the properties of SMGs 
with other g alaxie s at simi lar redshifts as prese nted by 
iDaddi et al.l (120071. Fi g 17) |Kajisawa_gt^alJ (|20ld . Fig. 2) 
and iRodighiero et all (120111. Fig . 1). For the SMGs we 
adopt the SFRs derived by MoT from the mid-far infrared 
emission. The comparison is shown in Figure [5J Here, for 
the nine different determinations of SMG masses considered 
in this paper, we plot SFR versus M* for all the individual 
SMGs in the sample, and compare the results with the loca- 
tion of the s o-called 'main-sequence' of star-forming ga lax- 
ies at z ~ 2 (|Daddi et alJl2007t IRodighiero et aLll201lft . 

In order to ensure that this comparison is fair and self- 
consistent, we have also investigated how the choice of SFH 
influences the location of this 'main-sequence'. We assumed 
that the SFR estimate is independent of the adopted SFH 
(as it is usually derived from bolometric or monochro- 
matic luminosity), and derived the s tellar masses of all 
galaxies the UPS field in CANDELS (jGrogin et al.l l201lt 
iKoekemoer et al.1 l201lh at z = 1.5-3 and M* = 10 9 5 - 
10 11 - 6 M (Cirasuolo et al. in prep.) using the same SED fit- 
ting method outlined in Section [2^21 for double-burst, single- 
burst and tau SFHs. If we adopt the results of the single- 
burst models for these galaxies, we obtained masses consis- 
tently lower by 0.1 dex relative to those derived using the 
tau models; on Fig. [S] (T burst' panels) this is indicated by 
the position of the coloured thick line towards lower stellar 
masses. On the other hand, as expected, the adoption of the 
double-burst models yields higher stellar masses. However, 
this shift transpires to be mass-dependent; at 10 10 M the 
double-burst SFH produces stellar masses 0.3 dex higher 
than the tau SFH, whereas this shift is only 0.1 dex at 
10 11 - 5 M (resulting in steepening of the thick line in the 
'double' panels). The latter is a significantly smaller shift 
than the 0.3 dex difference found between the double-burst 
and tau model stellar masses for SMGs (Tab. [IJ. 

It is clear from this figure that, for almost all models 
the SMGs lie essentially exactly on this 'main-sequence', 
within the errors. It is also clear that, even for the lower 
mass determinations, the average deviation from the 'main- 
sequence' is relatively modest. As expected (given Table [T] 
and the discussion in Section 2.2), the offset is g reatest for 
the ex ponentially-decreasing tau SFH with the iMarastonl 
(12005ft models, but even then it can be seen that the ma- 
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• M10 (cont+burst) 
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□ C10 (double) 
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C10 (1 burst) 
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Fig. 5. SFRs of SMGs as a fun ction of stellar m asses derived using the stellar population synthesis models of 
Bruzual fc Chariot] (|2003l left) and iMarastonl (|2005l . right). The symbols correspond to the ste llar masses derived by 
(red circles) . Irlill (blue stars) and our new estimates based on the SED modelling method of ICirasuolo et al"] i|2010l 



ClOf ) for three alternative forms of star formation history (SFH): double-component model (red squares), exponentially- 



declining (tau) model (light blue triangles) and single burst (green pluses). The median values are shown as large black 



Daddi et al 



Daddietal 



symb ols with 3a erro r bars. The 'main-seque nce' of star-forming galaxies at z ~ 2 is shown as a solid line 
l2007l) and dotted line (|Rodighiero et aT]l201l|) . The thick, coloured solid lines represent the 'main-sequence' of 
(|2007f ) at z ~ 2 corrected for the systematic difference in average stellar mass introduced by the self-co nsistent adoption 
of SFH to model SMGs in each panel (see Section 2] for details). A ll results have be en converted to the IChabri cr (2003) 
IMF. Unless the exponentially declining (tau model) SFH with the IMarastonl (|2005T ) models are assumed, the SMGs are 
completely consistent with the high- mass end of this sequence. This agreement is essentially perf ect for the Midi an d 
("rouble -component models. The magenta star shows the position of an average SMG as plotted bv lDaddi et all (|2007l) . 
Clearly, regardless of the assumed evolutionary stellar model or SFH, this supposed typical value of sSFR for SMGs is 
in fact unrepresentative of real SMGs. 
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jority of the SMGs are still actually consistent with the 
'main-sequence', and the formal ~ 3c offset in average 
sSFR is in fact driven by a subset of extreme outliers. It 
can also be seen that none of the various determinations 
of SMG stellar masses yields a sSFR for a typica l SMG 
as extreme as that claimed by iDaddi et al.l (|2007t ). This 
discrepancy is likely because their stellar mass estimates 
are based on dy namical masses from high-excitation CO 
line observations ()Tacconi et al.|[2006( ). which may bias the 
stellar mass estimates low, as we explained in Section [XT] 
Moreover, their use of infrared luminosity of ~ 10 13 L^, 
estimated from a single 850 (im detection (|Greve et al.l 
120051: iTacconi et al.l [2006h . represents an overestimate, as 
revealed by recent He rschel studies (|Magnelli et al.l [2012t 
iRoseboom et~aD I2TJI2I ~ lO 1 ^- 5 - 12 - 7 !^ making the SFR 
estimate of IDaddi et al.l (120071 ) a factor of 2-3 too high. 

To quantify the apparent consistency between the prop- 
erties of SMGs and other star-forming galaxies at compa- 
rable epochs we report, in Table [TJ the mean and median 
sSFRs for the SMGs as derived from the various alternative 
models (as stated above these values are I MF- indep endent ) . 
Again, with the possible exception of the IMarastonl (|2005l ) 
tau-model, all the median values ar e clearly consistent with 
the values of a few Gy r -1 found bv lDaddi et alj (|2007f ) and 
iKaiisawa et al.l (|2010f ). Hence, it would appear that, while 
a few SMGs may be caught at unusually spectacular mo- 
ments, in general it seems hard to avoid the conclusion that 
the SMG population is perfectly consistent with the high- 
mass end of the general star-forming galaxy population at 
these redshifts. Interestingly, it is only really the double- 
component models which yield mean and median values 
which arc formally consistent with each other, with the 
other models always producing a few extreme outliers with 
very low stellar masses. Thus, the use of the most physi- 
cally plausible star-formation histories yields consistent and 
fairly well defined typical mean sSFR values for SMGs of 
2 — 3Gyr _1 , in rather good agreement wit h the height of 
the pl ateau in sSFR at z > 2 discussed by iGonzalez et al.1 

5. Conclusions 

We have explored the potential origins of systematic dif- 
ferences in the stellar masses derived for SMGs, and have 
fo und t hat the d ifference between the stellar-mass estimates 
of lMld andlHTll can be explained by the combination of dif- 
ferent assumed IMFs, stellar population models and types 
of star formation histories. The available evidence supports 
the use of the Chabrier IMF (with the Salpeter IMF yield- 
ing excessive masses for several SMGs). The choice of which 
stellar synthesis model to adopt turns out to be relatively 
unimportant compared to the choice of star formation his- 
tory. The one exception to this is the combination of an ex- 
pone ntially decl i ning s tar formation history with the mod- 
els of IMarastonl (|2005l i. which produces radically different, 
unrealistically low stellar masses, presumably due to the 
very strong assumed TP-AGB contribution in these mod- 
els. 

For the same IMF and evolutionary synthesis model, 
the remainin g diff erenc e bet ween the average stellar mass 
estimates of iMlOl and iHllI can be accounted for by dif- 
ferences in adopted star-f orma tion history. Specifically, the 
two-component models of lMldl yield values ~ 3 tim es larger 
than the single-component models utilised bv lHlll . We have 



explored this issue in more detail, by running further single- 
burst, tau-mo del a nd double-compone nt model fits using 
the m etho d of ClOl and usi ng both the iBruzual &: Chariot! 



(120031 ) and IMarastonl (|2005f) models. We find that the mod- 
els which allow past and "present" star-formation activ- 
ity to be disconnected yield very consistent average stel- 
lar masses in the range M* = 2-3 x 10 11 M Q , regardless 
of precisely how this is parameterized. We arg ue, th ere- 
fore, that the main r eason that authors such as IHllI and 
IWardlow et al.l (|201ll ) have concluded in favour of typical 
SMG stellar masses M* < 1 x 10 n M© is that they have 
assumed models in which the SFR in SMGs was always as 
high or higher prior to the epoch at which the object has 
been selected on the basis of its extreme star-formation ac- 
tivity. The more physically realistic two-component models 
are supported by other evidence for rising star-formation 
histories in star- forming galaxies at z > 2. 

Considering potential orientation effects, and updated 
estimates of dynamical ma sses, w e find that even the higher 
stellar mass estimates of IMlOl with the Chabrier IMF 
certainly cannot be excluded with current CO line data, 
and indeed yield total masses in good agreement with dy- 
namical mass estimates. We also argue that SMG stellar 
masses M* ~ 2 x 10 11 M appear to be perfectly plausible 
given current measurements of the evolving mass function 
of the general galaxy population selected at near-infrared 
and mid-infrared wavelengths. 

We also show that the mass estimates of lMlOl are in good 
agreement with the predictions of hydrodynamical merger 
m odels of SMGs and the cold inflow scenario, whereas those 
of IH1 ll are closer to (but still not as low as) the much lower- 
mass predictions of certain semi-analytic models of galaxy 
formation. 

Finally, we find that the inferred specific star-formation 
rates of SMGs are perfectly consistent with those derived 
for other star-forming galaxies at comparable redshifts, sug- 
gesting that most SMGs simply represent the high-mass ex- 
tension of the 'main-sequence' of star-forming galaxies at 
z > 2, and in general should not be considered as outliers 
in the stellar mass-SFR relation. 
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Table A.l. Stellar masses of individual SMGs calculated using various methods, stellar population models and star 
formation histories. 



M, (logM )~ 



SMG 



z 



MKT 
Padova 
double 



Hll" 
BC03 
single 



Hll c 
M05 
single 



CIO 
BC03 
lburst 



C10 e 
BC03 
tail 



CIO 7 
BC03 
2burst 



CIO 9 
M05 
lburst 



CIO" CIO* 
M05 M05 
tau 2burst 



SMMJ030226. 17+000624.5 0.080 10.275 -99.000 -99.000 10.010 10.000 10.330 9.810 9.810 10.590 

SMMJ030227. 73+000653. 5 1.408 10.915 11.180 11.009 11.180 11.190 11.750 11.320 11.320 11.770 

SMMJ030231. 81+001031. 3 1.316 10.205 9.080 9.209 9.330 9.740 10.050 9.550 9.690 10.120 

SMMJ030236. 15+000817.1 2.435 11.285 10.880 10.689 11.180 11.170 11.310 11.250 11.200 11.570 

SMMJ030238. 62+001106. 3 0.276 8.525 10.790 10.639 8.220 8.520 9.020 8.520 7.200 9.480 

SMMJ030244. 82+000632. 3 0.176 10.205 9.830 9.559 10.100 10.090 10.390 9.950 9.950 10.410 

SMMJ105151. 69+572636.0 1.147 11.175 11.450 11.109 11.370 11.370 11.370 11.520 11.520 12.140 

SMMJ105155. 47+572312. 7 2.686 11.235 10.580 10.419 10.880 10.460 11.530 11.180 8.650 11.570 

SMMJ105158. 02+571800. 2 2.239 11.675 11.480 11.259 11.320 11.340 11.320 11.600 11.260 12.020 

SMMJ105200. 22+572420. 2 0.689 10.945 9.870 9.759 9.870 9.440 9.870 10.170 7.630 11.140 

SMMJ105201. 25+572445. 7 2.148 11.115 10.470 10.349 11.020 11.040 11.020 11.020 10.220 11.540 

SMMJ105207. 49+571904.0 2.689 11.655 11.010 10.879 11.190 11.230 11.460 11.050 11.050 11.510 

SMMJ105225. 79+571906. 4 2.372 11.385 11.070 10.879 11.660 11.650 11.660 11.760 11.760 12.010 

SMMJ105227. 58+572512. 4 2.142 11.565 11.370 11.189 11.190 11.420 11.190 11.310 11.190 11.750 

SMMJ105227. 77+572218. 2 1.956 11.445 9.820 9.609 -99.000 -99.000 11.190 -99.000 -99.000 -99.000 

SMMJ105230. 73+572209. 5 2.611 11.875 11.470 11.289 11.320 10.890 12.130 11.690 9.040 11.970 

SMMJ105238. 19+571651.1 1.852 10.905 10.700 10.519 11.160 11.460 11.680 11.410 11.210 11.480 

SMMJ105238. 30+572435. 8 3.036 11.975 11.450 11.289 11.180 10.820 11.180 11.500 9.000 11.820 

SMMJ123549. 44+621536. 8 2.203 11.425 11.320 11.149 10.820 10.830 10.820 10.880 10.880 11.520 

SMMJ123553. 26+621337. 7 2.098 11.075 10.650 10.489 10.640 11.290 11.500 10.950 10.770 11.310 

SMMJ123555. 14+620901. 7 1.875 11.505 11.230 11.069 11.260 11.740 12.090 11.550 11.150 11.900 

SMMJ123600. 10+620253.5 2.710 11.285 -99.000 -99.000 9.440 9.300 9.440 9.720 9.000 12.150 

SMMJ123600. 15+621047.2 1.994 11.385 10.470 10.249 11.140 11.520 11.650 11.410 11.630 11.840 

SMMJ123606. 72+621550. 7 2.416 11.195 10.450 10.369 11.010 11.500 11.650 11.330 10.810 11.720 

SMMJ123606. 85+621021. 4 2.509 11.735 11.270 11.149 11.370 11.500 11.720 11.340 11.340 11.750 

SMMJ123616. 15+621513. 7 2.578 11.715 10.760 10.799 11.280 10.860 11.280 11.580 9.040 12.050 

SMMJ123618. 33+621550. 5 1.865 11.445 10.880 10.659 10.770 11.340 11.630 11.090 10.900 11.430 

SMMJ123621. 27+621708. 4 1.988 11.655 10.780 10.619 10.890 11.460 11.760 11.210 10.710 11.570 

SMMJ123622. 65+621629. 7 2.466 11.665 11.250 11.059 11.020 11.600 11.890 11.340 11.210 11.720 

SMMJ123629. 13+621045.8 1.013 11.445 11.030 10.909 10.870 11.150 11.500 11.140 10.810 11.530 

SMMJ123632. 61+620800.1 1.993 11.075 10.520 10.359 11.060 10.670 11.060 11.370 8.830 11.740 

SMMJ123634. 51+621241.0 1.219 11.405 10.940 10.769 10.850 10.500 10.850 11.180 8.610 11.520 

SMMJ123635. 59+621424.1 2.005 11.125 11.420 11.239 11.610 12.270 12.470 11.920 11.600 12.270 

SMMJ123636. 75+621156.1 0.557 10.215 9.930 9.779 9.720 10.030 10.360 10.040 8.480 11.040 

SMMJ123651. 76+621221. 3 0.298 9.875 -99.000 -99.000 9.430 9.590 9.940 9.510 8.650 10.230 

SMMJ123707. 21+621408.1 2.484 11.695 11.260 11.139 11.050 11.090 11.050 11.370 10.570 11.750 

SMMJ123711. 98+621325. 7 1.992 11.035 10.610 10.419 10.560 11.150 11.420 10.880 10.560 11.230 

SMMJ123712. 05+621212. 3 2.914 11.845 11.340 11.139 11.030 11.660 11.800 11.350 11.480 11.680 

SMMJ123716. 01+620323. 3 2.037 11.675 -99.000 -99.000 12.010 11.850 13.150 12.100 10.000 13.100 

SMMJ123721. 87+621035. 3 0.979 11.185 10.800 10.689 10.720 11.130 11.530 10.910 10.780 11.780 

SMMJ131201. 17+424208.1 3.405 11.105 11.010 10.909 11.190 10.860 11.830 11.460 10.310 11.870 

SMMJ131208. 82+424129.1 1.544 11.245 10.980 10.809 10.770 10.470 10.770 11.060 9.800 11.440 

SMMJ131212. 69+424422. 5 2.805 11.565 10.710 10.649 10.970 10.960 10.970 10.830 10.830 11.420 

SMMJ131215. 27+423900. 9 2.565 12.065 11.980 11.869 12.530 12.100 12.530 12.820 10.290 13.220 

SMMJ131222. 35+423814.1 2.565 11.505 11.490 11.309 12.120 11.620 12.120 12.470 9.790 12.760 

SMMJ131225. 20+424344. 5 1.038 10.795 10.810 10.609 10.590 10.170 10.590 10.910 8.340 11.770 

SMMJ131225. 73+423941. 4 1.554 10.965 10.580 10.469 10.410 11.180 11.340 10.690 11.000 11.550 

SMMJ131228. 30+424454. 8 2.931 11.415 11.040 10.819 11.080 11.230 11.730 11.370 10.610 11.770 

SMMJ131231. 07+424609.0 2.713 10.995 -99.000 -99.000 15.890 15.330 13.360 12.320 12.320 12.450 

SMMJ131232. 31+423949. 5 2.320 11.765 10.180 10.089 12.180 12.450 12.800 12.460 12.200 12.720 

SMMJ131239. 14+424155. 7 2.242 11.255 10.700 10.529 10.530 11.020 11.140 10.780 10.380 11.340 

SMMJ141741. 81+522823.0 1.150 11.755 11.770 11.619 11.890 11.460 11.890 12.200 10.990 12.560 

SMMJ141742. 04+523025. 7 0.661 10.865 10.750 10.569 10.660 10.230 10.660 10.960 8.410 11.420 

SMMJ141750. 50+523101.0 2.128 10.845 10.010 9.829 9.900 10.360 10.940 10.010 9.610 11.040 

SMMJ141800. 40+512820. 3 1.913 12.185 11.160 10.999 12.450 13.620 12.450 12.740 11.790 12.740 

SMMJ141802. 87+523011.1 2.127 10.725 10.260 10.069 10.260 10.460 10.810 10.440 10.440 10.960 

SMMJ141809. 00+522803. 8 2.712 11.975 11.770 11.479 11.370 11.930 12.050 11.540 12.090 12.140 

SMMJ141813. 54+522923. 4 3.484 11.055 10.360 10.499 10.440 11.100 11.110 10.460 10.850 11.160 

SMMJ163627. 94+405811. 2 3.180 10.985 10.580 10.249 10.420 10.760 11.090 10.590 10.730 11.120 

SMMJ163631. 47+405546. 9 2.283 11.435 10.830 10.699 10.860 10.420 10.860 11.180 8.610 11.660 

SMMJ163639. 01+405635. 9 1.495 10.975 10.700 10.539 11.210 10.330 11.210 11.220 8.510 11.520 

SMMJ163650. 43+405734. 5 2.378 11.315 11.160 10.939 11.580 11.590 11.580 11.720 11.720 12.180 

SMMJ163658. 19+410523.8 2.454 11.515 11.110 10.959 11.050 11.750 11.910 11.350 11.700 11.760 

SMMJ163658. 78+405728.1 1.190 11.235 10.990 10.839 10.800 10.380 10.800 11.110 8.960 11.480 

SMMJ163704. 34+410530. 3 0.840 10.595 9.940 9.759 9.980 9.950 10.540 10.210 9.410 10.620 

SMMJ163706. 51+405313. 8 2.374 11.655 11.210 11.079 11.460 11.230 12.010 11.670 10.530 12.090 

SMMJ221724. 69+001242.1 0.510 11.055 -99.000 -99.000 10.580 10.150 11.680 10.890 8.340 11.900 

SMMJ221725. 97+001238. 9 3.094 11.705 -99.000 -99.000 10.860 11.480 11.470 11.170 11.380 11.510 

SMMJ221733. 02+000906.0 0.926 11.385 -99.000 -99.000 12.000 12.000 12.000 12.100 12.090 11.760 

SMMJ221733. 12+001120.2 0.652 11.075 10.700 10.579 10.420 10.550 10.420 10.740 10.140 11.680 

SMM.J221733. 91+001352.1 2.555 11.715 11.160 10.959 11.200 11.520 11.860 11.490 11.090 11.890 

SMMJ221735. 15+001537.2 3.098 11.115 10.710 10.459 10.620 11.050 11.280 10.900 10.970 11.280 

SMMJ221735. 84+001558. 9 3.089 11.155 10.860 10.689 10.930 11.350 11.720 11.220 11.070 11.600 

SMMJ221737. 39+001025.1 2.614 11.315 -99.000 -99.000 10.960 11.470 11.370 11.010 11.410 11.490 

SMMJ221804. 42+002154. 4 2.517 11.035 10.920 10.739 10.800 10.450 11.470 11.080 10.230 11.500 

SMMJ221806. 77+001245. 7 3.623 11.315 11.010 10.839 11.330 11.730 11.990 11.600 11.570 12.140 



Notes. The IChabrierl (|2003T ) IMF is assumed. -99 indicates the lack of an estimate for a given SMG. (a) Stellar mass derived 
bv iMld . (6) Stellar mass derived bv lHllI with BC03 models. (c) Stellar mass derived bv lHllI with M05 models. (d) Stellar mass 
derived as in CIO with BC03 models and sing le-burst SFH. (c) Stellar mass derived as in CIO with BC03 models and tau SFH. 

Stellar mass derived as in CIO with BC03 models and double-burst SFH. (s) Stellar mass derived as in CIO with M05 models 
d&d sing le-burst SFH. (h) Stellar mass derived as in CIO with M05 models and tau SFH. (i) Stellar mass derived as in CIO with 
M05 models and double-burst SFH. Re ference s: IM Id dMichalowski et al.ll2010af ): MTU flHainline et alj|2011f >: IClOl l|Cirasuolo et all 
[20101 ): BC03 (|Bruzual fc Charlot|[200l ): M05 (|Marastonll2005l ). 



Michalowski et al.: The stellar masses and specific star- formation rates of submillimetre galaxies 



Bouche N., Lehnert M.D., Peroux C, 2005, MNRAS, 364, 319 
Bressan A., Fagotto F., Bertelli G., Chiosi C, 1993, A&AS, 100, 647 
Bruzual G., Chariot S., 2003, MNRAS, 344, 1000 
Bussmann R.S., Dcy A., Armus L., et al., 2012, ApJ, 744, 150 
Cappellari M., Bacon R., Bureau M., et al., 2006, MNRAS, 366, 1126 
Caputi K.I., Cirasuolo M., Dunlop J.S., et al., 2011, MNRAS, 413, 
162 

Chabrier G., 2003, ApJ, 586, L133 

Chapman S.C., Blain A.W., Smail I., Ivison R.J., 2005, ApJ, 622, 772 
Cimatti A., Cassata P., Pozzetti L., et al., 2008, A&A, 482, 21 
Cirasuolo M., McLurc R.J., Dunlop J.S., et al., 2007, MNRAS, 380, 
585 

Cirasuolo M., McLure R.J., Dunlop J.S., et al., 2010, MNRAS, 401, 
1166 

Conroy C, Gunn J.E., 2010, ApJ, 712, 833 

Coppin K., Chapin E.L., Mortier A.M. J., et al., 2006, MNRAS, 372, 
1621 

Daddi E., Dickinson M., Morrison G., et al., 2007, ApJ, 670, 156 
Dave R., Finlator K., Oppenheimer B.D., et al., 2010, MNRAS, 404, 
1355 

Downes D., Solomon P.M., 2003, ApJ, 582, 37 
Drory N., Bender R., Feulner G., et al., 2004, ApJ, 608, 742 
Dunlop J.S., McLure R.J., Yamada T., et al., 2004, MNRAS, 350, 769 
Dye S., Eales S.A., Aretxaga I., et al., 2008, MNRAS, 386, 1107 
Dye S., Eales S., Moncelsi L., Pascale E., 2010, MNRAS, 407, L69 
Eminian C, Kauffmann G., Chariot S., et al., 2008, MNRAS, 384, 
930 

Engel H., Tacconi L.J., Davies R.I., et al., 2010, ApJ, 724, 233 
Erb D.K., Steidel C.C., Shapley A.E., et al., 2006, ApJ, 646, 107 
Fagotto F., Bressan A., Bertelli G., Chiosi C, 1994a, A&AS, 104, 365 
Fagotto F., Bressan A., Bertelli G., Chiosi C, 1994b, A&AS, 105, 29 
Fagotto F., Bressan A., Bertelli G., Chiosi C, 1994c, A&AS, 105, 39 
Finlator K., Oppenheimer B.D., Dave R., 2011, MNRAS, 410, 1703 
Fontanot F., Monaco P., Silva L., Grazian A., 2007, MNRAS, 382, 
903 

Genzel R., Baker A.J., Tacconi L.J., et al., 2003, ApJ, 584, 633 
Gerhard O., Kronawitter A., Saglia R.P., Bender R., 2001, AJ, 121, 
1936 

Girardi L., Bressan A., Chiosi C, Bertelli G., Nasi E., 1996, A&AS, 
117, 113 

Girardi L., Bressan A., Bertelli G., Chiosi C, 2000, A&AS, 141, 371 
Gonzalez J.E., Lacey C.G., Baugh CM., Frenk C.S., 2011, MNRAS, 
413, 749 

Gonzalez V., Labbe I., Bouwens R.J., et al., 2010, ApJ, 713, 115 
Greve T.R., Bertoldi F., Smail I., et al., 2005, MNRAS, 359, 1165 
Grillo C, Gobat R., Lombardi M., Rosati P., 2009, A&A, 501, 461 
Grogin N.A., Kocevski D.D., Faber S.M., et al., 2011, ApJS, 197, 35 
Hainline L.J., 2008, Multi- Wavelength Properties of Submillimeter- 

Selected Galaxies, Ph.D. thesis, California Institute of Technology 
Hainline L.J., Blain A.W., Smail I., et al., 2009, ApJ, 699, 1610 
Hainline L.J., Blain A.W., Smail I., et al., 2011, ApJ, 740, 96 
Hatsukade B., Iono D., Akiyama T., et al., 2010, ApJ, 711, 974 
Hayward C.C., Keres D., Jonsson P., et al., 2011a, ApJ, 743, 159 
Hayward C.C., Narayanan D., Jonsson P., et al., 2011b, Astronomical 

Society of the Pacific Conference Scries, 440, 369 
Hayward C.C., Jonsson P., Keres D., et al., 2012, MNRAS, submitted, 

arXiv: 1203. 1318 

Henriques B., Maraston C, Monaco P., et al., 2011, MNRAS, 415, 
3571 

Hezaveh Y.D., Holder G.P., 2011, ApJ, 734, 52 

Iglesias-Paramo J., Buat V., Hernandez-Fernandez J., et al., 2007, 
ApJ, 670, 279 

Ikarashi S., Kohno K., Aguirre J.E., et al., 2011, MNRAS, 415, 3081 
Ivison R.J., Smail I., Le Borgne J.F., et al., 1998, MNRAS, 298, 583 
Ivison R.J., Papadopoulos P.P., Smail I., et al., 2011, MNRAS, 412, 
1913 

Kajisawa M., Ichikawa T., Tanaka I., et al., 2009, ApJ, 702, 1393 
Kajisawa M., Ichikawa T., Yamada T., et al., 2010, ApJ, 723, 129 
Kncib J., van der Werf P.P., Kraiberg Knudsen K., et al., 2004, 

MNRAS, 349, 1211 
Knudsen K.K., Kneib J., Richard J., Petitpas G., Egami E., 2010, 

ApJ, 709, 210 

Koekemoer A.M., Faber S.M., Ferguson H.C., et al., 2011, ApJS, 197, 
36 

Kriek M., Labbe I., Conroy C, et al., 2010, ApJ, 722, L64 
Kroupa P., 2001, MNRAS, 322, 231 

Lima M., Jain B., Devlin M., Aguirre J., 2010, ApJ, 717, L31 
MacArthur L.A., McDonald M., Courteau S., Jesus Gonzalez J., 2010, 



ApJ, 718, 768 

Magnclli B., Lutz D., Santini P., et al., 2012, A&A, 539, A155 
Maraston C, 2005, MNRAS, 362, 799 

Maraston C, Daddi E., Renzini A., et al., 2006, ApJ, 652, 85 
Maraston C, Pforr J., Renzini A., et al., 2010, MNRAS, 407, 830 
Menendez-Delmestre K., Blain A.W., Alexander D.M., et al., 2007, 
ApJ, 655, L65 

Menendez-Delmestre K., Blain A.W., Smail I., et al., 2009, ApJ, 699, 
667 

Michalowski M., Hjorth J., Watson D., 2010a, A&A, 514, A67 
Michalowski M.J., Watson D., Hjorth J., 2010b, ApJ, 712, 942 
Mosleh M., Williams R.J., Franx M., Kriek M., 2011, ApJ, 727, 5 
Motohara K., Takata T., Iwamuro F., et al., 2005, AJ, 129, 53 
Murphy E.J., Chary R.R., Alexander D.M., et al., 2009, ApJ, 698, 
1380 

Narayanan D., Cox T.J., Hayward C.C., Younger J.D., Hernquist L., 

2009, MNRAS, 400, 1919 
Narayanan D., Dey A., Hayward C.C., et al., 2010a, MNRAS, 407, 

1701 

Narayanan D., Hayward C.C., Cox T.J., et al., 2010b, MNRAS, 401, 
1613 

Negrello M., Hopwood R., De Zotti G., et al., 2010, Science, 330, 800 
Noeske K.G., Weiner B.J., Faber S.M., et al., 2007, ApJ, 660, L43 
Paciga G., Scott D., Chapin E.L., 2009, MNRAS, 395, 1153 
Papovich C, Finkelstein S.L., Ferguson H.C., Lotz J.M., Giavalisco 

M., 2011, MNRAS, 412, 1123 
Pope A., Chary R.R., Alexander D.M., et al., 2008, ApJ, 675, 1171 
Pozzetti L., Bolzonella M., Lamareille F., et al., 2007, A&A, 474, 443 
Ricciardelli E., Trujillo I., Buitrago F., Conselice C.J., 2010, MNRAS, 

406, 230 

Riechers DA., Hodge J., Walter F., Carilli C.L., Bertoldi F., 2011, 
ApJ 

Rodighiero G., Daddi E., Baronchelli I., et al., 2011, ApJ, 739, L40 
Roseboom I.G., Ivison R.J., Greve T.R., et al., 2012, MNRAS, 419, 
2758 

Salpeter E.E., 1955, ApJ, 121, 161 

Santini P., Maiolino R., Magnelli B., et al., 2010, A&A, 518, L154 
Schael A.M., 2009, The Star-Formation History of Massive Galaxies, 

Ph.D. thesis, University of Edinburgh 
Silva L., Granato G.L., Bressan A., Danese L., 1998, ApJ, 509, 103 
Swinbank A.M., Laccy C.G., Smail I., et al., 2008, MNRAS, 391, 420 
Swinbank A.M., Smail I., Chapman S.C., et al., 2010a, MNRAS, 405, 

234 

Swinbank A.M., Smail I., Longmore S., et al., 2010b, Nature, 464, 733 
Tacconi L.J., Neri R., Chapman S.C., et al., 2006, ApJ, 640, 228 
Tamura Y., Iono D., Wilner D.J., et al., 2010, ApJ, 724, 1270 
Targett T.A., Dunlop J.S., McLure R.J., et al., 2011, MNRAS, 412, 
295 

Taylor M.B., 2005, In: P. Shopbell, M. Britton, & R. Ebert (ed.) 
Astronomical Data Analysis Software and Systems XIV, vol. 347 
of Astronomical Society of the Pacific Conference Series, 29 
Thomson M.G., Chary R.R., 2011, ApJ, 731, 72 
Valiante E., Lutz D., Sturm E., et al., 2007, ApJ, 660, 1060 
van der Wei A., Franx M., van Dokkum P.G., et al., 2006a, ApJ, 636, 
L21 

van der Wei A., Franx M., Wuyts S., et al., 2006b, ApJ, 652, 97 
van Dokkum P.G., 2008, ApJ, 674, 29 
van Dokkum P.G., Conroy C, 2010, Nature, 468, 940 
van Dokkum P.G., Conroy C, 2011, ApJ, 735, L13 
Vieira J.D., Crawford T.M., Switzer E.R., et al., 2010, ApJ, 719, 763 
Wardlow J.L., Smail I., Coppin K.E.K., et al., 2011, MNRAS, 415, 
1479 

Watabc Y., Risaliti G., Salvati M., et al., 2009, MNRAS, 396, LI 
Weidner C, Kroupa P., Pflamm-Altenburg J., 2011, MNRAS, 412, 
979 

Wilson G.W., Hughes D.H., Aretxaga I., et al., 2008, MNRAS, 390, 
1061 

Zibetti S., Chariot S., Rix H., 2009, MNRAS, 400, 1181 
Zibetti S., Gallazzi A., Chariot S., Pasquali A., Pierini D., 2012, 
arXiv: 1203.4571 



13 



